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Abstract 

We establish an analog of the Beale-Kato-Majda criterion for singularities of smooth 
solutions of the system of PDE arising in the Oldroyd model for ideal viscoelastic flow. 

It is well known that smooth solutions of the initial value problem associated with Euler's 
equations 

^ ^ ^ , {x,t)eR'x{0,T) [1.1] 
V ■ u = 

exist for some finite time T > 0. Here u = u{x,t) G M.^ and p = p{x,t) G M represent 
the local velocity and pressure of a perfect fluid, respectively. One of the major challenges 
in PDE theory is to deduce whether or not finite time singularities do indeed occur A 
celebrated result of Beale, Kato, and Majda asserts that if a solution u of Euler's equations 
possesses a singularity at a finite time T, then necessarily 

|(V X u)(t)\L^(M.3)dt = +00 



[1]. In particular, the L°°(]R3) norm of the vorticity must blow up as t approaches T from 
below 

limsup |(V X M)(t)|ioo(]R3) = +00. 

Along with relatively standard energy estimates for u and V x u, the key insight for 
Beale, Kato and Majda was a clever application of Kato's inequality: 

|V/|loo(k3) < C {1 + (1 + ln+ |/U3(M3))|V X /|l.o(u3) + |V X /|l2(k3)} (1.2) 



*This material is based upon work supported by the National Science Foundation under Grant No. DMS- 
1004733. 

singularity is a time T > where u{T) fails belong to a natural function space such as C°°(M'^;R^). 



1 



for all / e H^{M.^;R^) with 

V-/ = 



[3]. Here C is a universal constant, and 

In^ X = 



Inx, X > 1 
0, X < 1 



The above inequality was derived from a careful analysis of the Biot-Savart law 

fix) = I V0(x - y) X (V X f{y))dy, x G M^ 



where (j){x) = l/47r\x\ is the fundamental solution of Laplace's equation on M^. 

Caflisch, Klapper, and Steel used these ideas to prove a similar blow up criterion for the 
ideal MHD equations 

'dtu + {u ■ V)n = -V{p+ \b\y2) + {b ■ V)b 
dtb+{u-V)b= {b-V)u , (x, t) G x (0, T). (1.3) 

V-M = V-6 = 

They showed that if the solution to the inital value problem, with smooth divergence free 
data, has a singularity at time T, then it must be that 

T 

|(V X n)(t)|ioo(K3) + |(V X 6)(t)|Loo(ig3)dt = +00 

JO 

0. 

In this short note, we develop these ideas further and establish an analogous blowup 
criterion for solutions of the following system of PDE 

dtU + {u ■ V)u = -Vp + V ■ FF* 
dtF + {u-V)F = VuF , (a;,t) G X (0,T). (1.4) 

V ■ n = 

This system of PDE arises in the Oldroyd model for ideal viscoelastic flow i.e. a viscoelastic 
fluid whose elastic properties dominate its behavior. Here F = F(x, t) G M^^^ represents the 
local deformation gradient of the fluid. See [7] for more on this model. 

We remark that global existence for solutions near equilibrium of the viscous analog of 
(11. 4p has been verified by Lin et.al. [7] for two-dimensional flow and Lei et.al. [5] for three- 
dimensional flow; see also the work Lin and Zhang P| for a further discussion of these topics. 
In these references, three assumed properties of the deformation gradient are crucial for the 
methods employed: 

• The condition V ■ F* = (see equation (II. 5p below); 
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The curl of the deformation gradient is of "higher order" (see equation (3) in [5]); 



• The condition detF = 1. 

Note that the third condition is natural as the flow is assumed to be incompressible. For 
the vanishing viscosity case, Sideris and Thomases [10] also established the global existence 
of smooth solutions near the equilibrium using the incompressible limit by imposing a null 
condition on the elastic stress. 

One motivation of this work is to consider which conditions in the above list are nec- 
essary to assume when considering global smooth solutions without any type of smallness 
requirement. The central result of this work suggests that even if the incompressibility of 
the deformation gradient is imposed, smooth solutions of (11. 4p may become singular in finite 
time. From this standpoint, the property that the curl of the deformation gradient is higher 
order is perhaps the most important when considering global existence. 



Before stating our main result, let us make a simple observation alluded to above. Letting 

Fk = Fck denote the columns of F we take the divergence the second equation in (ll.4p to 
arrive at 

9t(V ■ Fk) + (n ■ V)(V ■ Ffc) = 0, = 1, 2, 3. (1.5) 

Therefore, if V ■ -Ffc = initially, it will remain so for later times. In what follows, we will 
make this assumption and also note that it implies the equality 



V-FF' = J2i^k-V)Fk. 



k=l 

With this reduction, equation (11. 4p becomes 

dtU + {u ■ V)u = -Vp + ELi(^fc ■ V)Ffc 
dtFk + iu-\/)Fk = iFk-V)u (1.6) 

V ■ M = V ■ Ffc = 

for k = 1,2, 3. We are now ready to state and prove our main result. 
Theorem 1.1. Assume that 

uo e H%R^) and Fk,o G //'(M^) 
for k = 1, 2, 3 and some s G N, s > 3. Furthermore, assume 

V ■ Mo = V ■ Ffc,o = 0. 

Then equation (\1.6\\ has a solution u, Fk with respective initial data uq, F^a belonging to the 
space 

C([0, T], H'{R^)) n C\[0, T], H'~\R^)), (1.7) 
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provided 

/ |(V X n)(t)U^(K3) + V|(V X Ffc)(t)|L^(M3)rft < OO. (1.8) 

■^0 k=l 

In particular, solutions with smooth initial data remain smooth on (0, T) provided (11. 8p holds. 

Proof. The short time existence of solutions u, with respective initial data uq, F^ q belong- 
ing to the space (11.71) for some T > follows from routine arguments; it is straightforward 
to adapt the arguments in chapter 3, section 2 of 0, for instance. Moreover, if u, Fk do not 
belong to the space (11.71) then it must be that 

limSUp \u{t)\Hs(M.3) + ^ \Fkit)\H^(^E.S) = +00. 



t-fT 



k=l 



If this is not the case, there is Cq > such that \u{t)\Ho(M.s) + Ylk=i \^k{t)\H'>{R3) < Co for 
t < T. Moreover, letting ti < T we can solve (11.61) on the interval [ti, ti+T^) for some Tq > 
only depending on Cq and, in particular, independent of ti. By choosing ti G (T — Tq,T), we 
have extended the original solution past the time T. However, in this scenario u, Fk already 
belonged to the space (11.71) by the local existence of solutions. 
Therefore, it suffices to show that if (11.81) holds then 



Ht)\ls + J2\Fk{t)\% <C, 0<t<T (1.9) 

k=l 

for some constant C depending only on T, the H^{M.^) norms of the initial data, and 

„T 3 

M := |(V X M)(t)|Loo(R3) + V |(V X Fk)it)\Lo.(R3)dt < OO. 

•^0 k=l 

As in [1] and [2], the proof is accomplished in three steps: obtaining energy estimates on 
solutions, obtaining estimates for the curl of solutions, and applying the crucial estimate 
(II. 2p . The proof is virtually the same for any integer s > 3, so for simplicity we assume that 
8 = 3. For ease of notation, we will also denote L~(M3), L'^{R^), and H^{R^) by L°°, L^, and 
H^. We shall also write Ylk Yfk=i- Finally, all integrals will be taken over the entire 
space R^. 

Step 1 (Energy estimates). We first show that there is a universal constant C such that 



\u 



(t)lH3 + J2 \Fk{t)\l3 < (^uo\l, + ^ \Fk,o\i?j exp (^C \Vu\loo + ^ |VF; 




for < t < T. Let a be a multi-index with |a| < 3 and set 

V := ay and Gk := d^Fk 
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for A: = 1, 2, 3. Note that 

^3 



dtv + {u ■ V)v = -V{dy) + ELi(^fc ■ V)Gfe - Q 
dtGk + {u ■ V)Gk = {Fk ■ V)v - Rk 



where 



Q = d-iiu ■ V)u) - (u . V)9> - miF, . V)F,) - (F, ■ V)d-F,] 

Ru = ■ V)Ffe) - (n ■ V)<9:Ffc - ■ V)n) - (F^ ■ V)a>], A; = 1, 2, 3 



As V ■ n = V ■ -Ffc = 0, we have 



r \v\ 

di ^ 



2 

■ dtV 



Likewise, 



= -y" Y^Gk-iFk-V)v- jvQ 

k 

<-[ V)t; + |^U2|g|i2. 

k 

= I Y,-Gk-Rk + Gk-{Fu-V)v 

k 

<Y,\Gk\LARkW + I Gk-{Fk-V)t 



Consequently, 



Recall the identity, 

mfa) - fd^gy < G{\f\H-A9\Lo^ + |V/Uoo|^|^.} 
for smooth f,g and |a| < 3 |T]; here G is an absolute constant. This identity implies 
\Q\l^ < C{|mU3|VmUoo + J2k iFkln-A'^FklL^} 

\Rk\L^ <G{\u\HA^Fk\L^ + \Fk\HA^u\L^}, k = 1,2,3 



Hence, 



dt I ~Y 



k 

Summing over all multi indices |«| < 3 gives 



for some universal constant C. Inequality fll.lOp now follows from Gronwall's inequality. 
Step 2 (L^ estimates). Set 

w:='Vxu and := V x F^, k = 1,2,3. 

We claim 



Wli2 + X^|r,(t)|i2 < ^koli2 + ^|rfc,oli2jexp||c^ kU- + Ir^U^rfsj (1.11) 



for < t < T and some universal constant C. With hypothesis (11.81) . the above inequality 
implies an bound on w{t) and rk(t) that is independent of t G [0,T]. On our way to 
proving this bound, we recall a simple fact: u,v E iJ^(M'^; M'^), then 

V X [(u ■ V)v] = a(Vn, Vv) + (u ■ V)(V x t;) 

where a : M^^^ x M^^^ — )■ is given by 

3 

a{X, Y) := Y^iX'ei) x (Fe,), X, F e M^^^ 

i=l 

Taking the curl of the first two equations in (11.61) gives 

dtW + iu-V)w = j:,iFk-V)n-S 
_ dtrk + (n ■ V)rfc = {Ft ■ V)w - Tk 

where 

' S = a{Vu, Vu) - J2k «(Vi^fc, VFfc ) 
Tk = a{Vu,VFk)-a{VFk,Vu), k = 1,2,3. 



6 



We compute 



d f \w\ 
dt 



w ■ dfiu 



■ V)n - s] 

k J 



Likewise, 



k k 



Consequently, 



k 



< \\/u\l^\w\l^\Vu\l^ + ^ |VFfc|Loo|M;|i2|VFfc|L2 

k 

+ ^ (|Vn|Loo|rfc|L2|VFfe|i2 + | VF^I^oo |rfc|i2 1 Vu|2,2 



as 



|Vt;|L2 < C|V X v\l2 

for f G H^, V ■ f = 0. Inequality (II. lip now follows from Gronwall's inequality. 

Step 3 (Applying Kato's inequality). Now we are in a position to complete the proof. 
Assume that (11. 8p holds so that (II. lip implies w,rk G L°°((0, T); L^). By Kato's inequality 
(II. 2p . we have 

|Vu|L^+^|VFfc|L- <c|l+ (^w\l^ + ^\rk\L^^ ^ln(|u|/^3 + e) +ln ^^[Ffcl^a + ej j 
By the energy estimate (ll.lOp . we also have 

X (^n{\u\H, + e) + In \Fk\Hs + j IdsV 



\u\ 
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Now setting 



y{t) := In (|n(t)U3 + e) + In V \Fk{t)\H^ + e , < t < T 



y{t) <B + C [ I 1 + ( \w{s)\lo. + J2 \n{s)\L^ ] y{s) \ ds. 




gives 



Here C is a universal constant and B depends only on the norms of the initial data. By 
Gronwall's inquality, we conclude that y{t) is bounded by a constant depending on M, T 
and the norms of the inital data. This bound establishes inequality fll.9p and completes 
the proof of the theorem. □ 

Finally, we remark that results closely related to Theorem 11.11 have been obtained for two 
dimensional flow |6] and for three dimensional flow in the so-called "creeping flow" regime 
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